Heavy quarks at RHIC and LHC within a partonic transport model 



Jan UphofP, Oliver Fochler 3 , Zhe Xu b a , Carsten Greiner" 

" Institut fiir Theoretische Physik, Johann Wolfgang Goethe-Universitdt Frankfurt, Max-von-Laue-Str. I, D-60438 Frankfurt am Main, Germany 
''Frankfurt Institute for Advanced Studies, Ruth-Moufang-Str. 1, D-60438 Frankfurt am Main, Germany 



Abstract 

Production and space-time evolution of heavy quarks in central and non-central heavy-ion collisions at RHIC and LHC 
are studied with the partonic transport model Boltzmann Approach of MultiParton Scatterings (BAMPS). In addition 
to the initially created heavy quarks in hard parton scatterings during nucleon-nucleon collisions, secondary heavy 
quark production in the quark-gluon plasma is investigated and the sensitivity on various parameters is estimated. In 
BAMPS heavy quarks scatter with particles of the medium via elastic collisions, whose cross section is calculated 
with the running coupling and a more precise implementation of Debye screening. In this framework, we compute 
the elliptic flow and nuclear modification factor of heavy quarks and compare it to the experimental data. 
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1. Introduction 

Several experimental observations indicate that a medium of free quarks and gluons - the quark-gluon plasma 
(QGP) - is produced in ultra-relativistic heavy-ion collisions Charm and bottom quarks are an ideal probe for 

the early stage of these collisions since they can only be created in initial hard parton scatterings of nucleon-nucleon 
interactions or in the beginning of the QGP phase, where the energy density is still large. After their production 
they interact with other particles of the medium and can, therefore, reveal important information about the properties 
of the QGP. Flavor conservation renders them as an unique probe since they are tagged by their flavor even after 
hadronization. 

The experimental data of the elliptic flow V2 and nuclear modification factor Raa of heavy quarks ^ 0, 0] show 
that the energy loss of charm and bottom quarks is comparable to that of light quarks. Whether this large energy loss 
is due to collisional or radiative interactions - or both (or even other effects) - is under investigation (see 15|] for a 
recent overview and comparison with data). 

After the introduction of the parton cascade BAMPS we will discuss the production of heavy quarks at RHIC and 
LHC. In Sec.|4]our results on the elliptic flow and nuclear modification factor at RHIC are discussed and compared to 
the experimental data. 

2. Parton cascade BAMPS 

For the simulation of the QGP we use the partonic transport model Boltzmann Approach of MultiParton Scatterings 
(BAMPS) [6, 7], which describes the full space-time evolution of the QGP by solving the Boltzmann equation, 



3i + ^j f ^^ C " +C < + -' (1) 

for on-shell partons and pQCD interactions. Details of the model, the implemented processes, and the employed cross 
sections can be found in [ 6 , 7j, [stl . 
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Figure 1; Number of charm quark pairs produced in a central Pb+Pb collision at LHC with ysjvw = 5.5 TeV (left) and tJsnn = 2.76 TeV 
(right) according to BAMPS. The initial parton distributions are obtained with PYTHIA and the mini-jet model (and the color glass condensate for 
ysAW = 5.5 TeV). In all cases the initial charm quarks are sampled with PYTHIA for better comparison. 



3. Heavy quark production at RHIC and LHC 

In heavy-ion collisions charm and bottom quarks are produced in hard parton scatterings of primary nucleon- 
nucleon collisions or in the QGP. To estimate the initial heavy quark yield, we use PYTHIA )0] and scale from proton- 
proton collisions to heavy-ion collisions with the number of binary collisions. Secondary heavy quark production in 
the QGP is simulated with BAMPS. For the initial gluon distributions, the mini-jet model, the color glass condensate 
model and also PYTHIA are employed. 

In Au+Au collisions at RHIC with -sJJnn = 200 GeV between 0.3 and 3.4 charm pairs are produced in the QGP, 
depending on the model of the initial gluon distribution, the charm mass and whether a K — 2 factor for higher order 
corrections of the cross section is employed [8]. This is only a small fraction of the initially produced charm quarks 
and can be neglected for the most probable scenarios. At LHC with the much larger initial energy density, secondary 
charm production is enhanced and not negligible. In Pb+Pb collisions with ysAW - 5.5 TeV between 11 and 55 
charm pairs are produced during the evolution of the QGP [8] (see left panel of Fig. [T}. These values are of the same 
order as the initial yield. As is shown in the right panel of Fig.Q] even in the 2010 run with -\[snn = 2.76 TeV between 
5 and 28 charm pairs are created during the QGP phase. 

Bottom production in the QGP, however is very small both at RHIC and LHC and can be safely neglected. As a 
consequence, all bottom quarks at these colliders are produced in initial hard parton scatterings. 

Further details on heavy quark production can be found in Ref. Hfl. 



4. Elliptic flow and nuclear modification factor of heavy quarks at RHIC 

The elliptic flow and the nuclear modification factor 

-1 P AZA\ p d 2 N AA /dp T dy 

2 = \7TI ' = N bm d 2 N pp /dp T dy (2) 

(p x and p y are the momenta in x and y direction in respect to the reaction plane) of heavy quarks at mid-rapidity are 
observables which are experimentally measurable and reflect the coupling of heavy quarks to the medium. A large 
elliptic flow and a small Raa indicate strong interactions with the medium and a sizeable energy loss. Experimental 
results reveal that both quantities are of the same order as the respective values for light particles 13l|4j,|5|]. 

The leading order perturbative QCD cross section with a constant coupling a s = 0.3 and the Debye mass for 
the t channel screening is too small to build up the large elliptic flow measured at RHIC [10]. However, if we take 
the running of the coupling into account and determine the screening mass from comparison to hard thermal loop 
calculations, we obtain a \'i and Raa, which are much closer to the data. 
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Figure 2: Elliptic flow V2 (left) and nuclear modification factor Raa (right) of heavy quarks and heavy flavor electrons with pseudo-rapidity 
\rj\ < 0.35 at the end of the QGP phase for Au+Au collisions at RHIC with an impact parameter of b = 8.2 fm. The cross section of gQ — » gQ is 
multiplied with the factor K = 4. For comparison, data of heavy flavor electrons | 5] is shown. 



The following calculations are done analogously to fill IT2L IT3I [Till . An effective running coupling is obtained 
from measurements of <? + e~ annihilation and non-strange hadronic decays of t leptons 115111 ill . Since the t channel 



of the gQ — > gQ cross section is divergent, it is screened with a mass proportional to the Debye mass mo\ 

1 1 



t t — k m 2 D 



(3) 



The Debye mass is calculated by the standard definition m 2 D - An(\ + A///6) a s (f) T 2 , but with the running coupling 
for consistency. The prefactor k in Eq.[3]is mostly set to 1 in the literature without a sophisticated reason. However, 
one can fix this factor by comparing the energy loss per unit length AE/dx of the born cross section with k to the 



energy loss within the hard thermal loop approach to k » 0.2 01 1111211 . 

This more accurate treatment increases the cross section of elastic gluon heavy quark scattering by about a factor of 
10. Fig.|2]shows the elliptic flow V2 and nuclear modification factor Raa for heavy quarks and for heavy flavor electrons 
as a function of the transverse momentum pj. To yield the same values for these variables as the experimental data, 
the leading order cross section of elastic collisions with the running coupling and improved Debye screening must 
still be multiplied by a K factor of 4. We assume that this artificial K factor stands for the contribution of radiative 
energy loss. However, it must be checked if these corrections have indeed a similar effect as a constant K factor of 4. 
Therefore, the calculation of the next-to-leading order cross section is planned for the near future and will complement 
2 <-> 3 interactions for gluons, which are already implemented in BAMPS Jg]. 

Especially for high pj the shape of the V2 curve of heavy quarks is different from the experimental data. The 
reason for this discrepancy is that, experimentally, due to confinement, not heavy quarks, but heavy flavor electrons 
are measured. The latter stem from the decay of D and B mesons, which in turn are produced during hadronization of 
the QGP by merging of the charm or bottom quark with a light quark. However, despite the hadronization and decay 
processes heavy flavor electrons still reveal information about heavy quarks. Essentially, the shape of their spectrum 
is the same as for heavy quarks, but shifted to lower pj due to the decay process. 

For the description of the hadronization process of charm (bottom) quarks to D (B) mesons, we use Peterson 
fragmentation [ 16]. The decay to heavy flavor electrons is carried out with PYTHIA. Fig.|2]shows that the theoretical 
curves for heavy flavor electrons are in good agreement with the experimental data for high pj. For low pj, however, 
Peterson fragmentation is not a good description of the hadronization process and another scheme like coalescence 
must be employed. In the coalescence picture, the light quarks of the D/B mesons contribute also to its elliptic flow 
or nuclear modification factor, which increases both. Studies on the V2 and Raa of gluons in BAMPS are presented in 
1 3 HE III & 
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5. Conclusions 



In this talk we presented results on heavy quark production as well as on the elliptic flow and nuclear modification 
factor of heavy quarks. Charm production in the QGP at RHIC is to a good approximation negligible, but at LHC a 
significant fraction of the total charm number is produced in the QGP - even at a center-of-mass energy of yfsmf = 
2.76 TeV. Bottom production in the medium is negligible at RHIC and LHC. The leading order cross section of heavy 
quark scatterings with particles from the medium is too small to explain the experimentally measured elliptic flow and 
nuclear modification factor. However, a more precise implementation of Debye screening and the explicit running of 
the coupling enhances the cross section and yields results for heavy flavor electrons, which are much closer to the 
data, although a K factor of 4 must be employed for a good agreement with the data. In the future we will study if this 
simple multiplication of the cross section with a K factor can indeed account for higher order contributions. 
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